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A quark action designed for highly anisotropic lattice simulations is discussed. The mass- 
dependence of the parameters in the action is studied and the results are presented. Applications 
of this action in studies of heavy quark quantities are described and results are presented from 
simulations at an anisotropy of six, for a range of quark masses from strange to bottom. 



I. INTRODUCTION 

The anisotropic lattice has proved an invaluable tool 
for simulations of a variety of physical quantities. The 
precision calculation of the glueball spectrum was an 
early application of the approach |l| and it was recognised 
that anisotropic actions may also be advantageous in 
heavy quark physics calculations 2] . Correlators of heavy 
particles such as glueballs and hadrons with a charm or 
bottom quark have a signal which decays rapidly. Monte 
Carlo estimates of these correlation functions can be 
noisy, making it difficult to resolve a plateau over a con- 
vincing range of lattice time steps. Increasing the number 
of timeslices for which the effective mass of a particle has 
reached a plateau solves this problem and also decreases 
the statistical error in the fitted mass. Since this value 
may be used as an input to determine many physical pa- 
rameters this decrease is very beneficial. 

Secondly, improved precision in effective mass fits 
means that momentum-dependent errors of O(ap) can be 
disentangled from other discretisation effects and larger 
particle momenta may be considered. This is particularly 
relevant for the determination of semileptonic decay form 
factors where the overlap of momentum regions accessi- 
ble to experiments and lattice calculations is currently 
very small. Typically, experiments have more events with 
daughter particle momentum at or above 1 GeV. This is 
also the region where large momentum-dependent errors 
are expected in lattice calculations. The form factors of 
decays like B — > n£u and B — » if* 7 are inputs to deter- 
minations of CKM parameters so that increased preci- 
sion in lattice calculations can lead to tighter constraints 
on the Standard Model. This has motivated a study of 
2+2 anisotropic lattices where the temporal and one spa- 
tial direction are made fine and all momentum is injected 
along this fine spatial axis. Details of the progress to date 
in this work are in Refs. and Q. The 2+2 formula- 
tion has also proved useful for a precision determination 
of the static interquark potential over large separations, 
which is described in Ref. 0. In this paper we consider 
a 3+1 anisotropic action. The temporal lattice spacing, 
at is made fine relative to the spatial spacing, a s . The 
action is designed with simulations at large anisotropics 
in mind. To simulate a bottom quark with a relativis- 
tic action requires a lattice spacing of less than 0.04 
fm which is prohibitively expensive on an isotropic lat- 



tice where the simulation cost scales at least as (D(a~ A ). 
The anisotropic lattice offers the possibility of relativistic 
heavy-quark physics using reasonably modest computing 
resources. In the rest frame of a hadron with a heavy 
constituent, the quark four- momentum is closely aligned 
with the temporal axis, allowing an anisotropic discreti- 
sation to represent accurately the Dirac operator on the 
quark field. 

Implementing an anisotropic programme incurs a num- 
ber of computational overheads not associated with the 
isotropic lattice. The ratio of scales, £ = a s /a t deter- 
mined by studying a physical long-distance probe de- 
pends on bare parameters in the lattice action. While 
this dependence is straightforward to establish at the 
tree-level of perturbation theory, quantum corrections 
can occur at higher orders. In the quenched approxima- 
tion on a 3+1 lattice this is not a serious additional cost 
as the tuning can be done post-hoc. More worringly, in 
Ref. |5j it was pointed out that the choice of £ and its re- 
lation to the Wilson parameter, r on anisotropic versions 
of the Sheikholeslami-Wohlert (SW) action could intro- 
duce 0{a s m q ) errors. It is exactly errors of this type that 
the anisotropic action seeks to avoid and the appearance 
of these terms represents a serious tuning problem. 

In this paper we use an action, specifically designed 
for highly anisotropic lattices ie. £ > 5. By applying 
different improvement terms in the spatial and temporal 
directions the action is both doubler and 0(a s m q ) error 
free. This opens up the possibility of simulating directly 
at the bottom quark mass using a relativistic action. In 
addition, in this feasibility study the speed of light was 
determined at « 1% accuracy. This precision was gov- 
erned by finite statistics and could certainly be improved 
upon. 

The paper is organised as follows. The construction of 
the action is described in Section^] Scction lTTII compares 
this action with the sD34 action proposed in Ref. [g and 
details some analytic results. Results from a study of 
the dispersion relations and the mass-dependence of the 
speed of light are described in Section ITVl Some prelimi- 
nary results of this study have appeared in Ref. |jj . Our 
conclusions and a discussion of future work are contained 
111 Section |V| 
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II. DESIGNING HIGHLY ANISOTROPIC 
ACTIONS 

We begin by considering a Wilson-type action with 
Symanzik improvement to remove discretisation errors. 
Full O(a)-improvement requires a clover term and a field 
rotation, given by 

TO 

i, = [l--qp- m )W, (l) 
i, = ^[l- r ±gp- m )l (2) 

where a is the lattice spacing on an isotropic lattice and 
r is the usual Wilson parameter. The rotation described 
by Eqs. Q and (J2J) preserves locality and maintains a 
positive transfer matrix so that ghost states do not arise 
in a calculation of the free fermion propagator. However, 
in an anisotropic implementation of this action, when 
a t is made very small, these rotations may lead to the 
reappearance of doublers, an undesirable side-effect of 
the anisotropy. 

We would like to maintain the useful properties of ac- 
tions with nearest neighbour temporal interactions only. 
In particular, the positivity of the transfer matrix guar- 
antees that effective masses approach a plateau from 
above. Therefore, to construct an action suitable for high 
anisotropies we begin by applying field rotations in the 
temporal direction only, rewriting Eqs. and @ as 

V> = [l-^^oAj-mM, (3) 
$ = V?[l-^( 7 oA)-m)]. (4) 

This leads to a new action in which the temporal and 
spatial directions are treated differently. Having applied 
the rotations of Eqn. @ the continuum action is given 
by 

S' = $M r $ - (Dl - 9 - ei E t ) (5) 

where M r = /J, r JiDi+joDo + and /i r = (1 + ^rajm). 
At the tree-level, the rotations described in Eqs. © and 
Q do not generate a spatial clover term. As a result the 
(a ■ B) term does not appear in Eq. J2J). The chromo- 
electric field, Ei is defined as 

igE i = [D i ,D ], (6) 

and ei = <t m is given by ^ = — [ji, Jo]- 

ii 

The temporal doublers are removed by discretising the 
Dq term in the usual way. However, with no spatial ro- 
tation the spatial doublers remain and must be treated 
separately. They are removed by adding a higher-order, 
irrelevant operator to the action. This was first suggested 
by Hamber and Wu in Ref. The simplest such oper- 
ator is a spatial D 4 term which is added ad hoc to the 



Dirac operator giving an action, 

S' = ^M^ - (Dl - ^ciEi) if + saty £ Dty. 

i 

(7) 

This approach has previously been discussed in detail 
in Ref. In this formulation, s is a Wilson-like pa- 
rameter which is chosen such that the doublers receive 
a sufficiently large mass. The discretisation of the ac- 
tion in Eq. J7J is now straightforward. Only the jiDi 
term requires an improved discretisation since the sim- 
plest discretisation would lead to 0(a 2 ) errors. For this 
case we write 

= l{~[<Kx + a)-<Kx-a)] 

-L[t(x + 2a)-<f>(x-2a)]\ (8) 

and similarly the (unimproved) discretisations of <9, d 2 
and <9 4 are 

A« = J_{^( x + o)-^(i-o)}, (9) 
la 

bS^Mx) = \ + a) + (p(x - a) - 20} , (10) 
a 

A( 4 )0(x) = \{[<f>(x + 2a) + (f>(x-2a)} 
a 1 

- 4[(/>(x + a)-<j>(x-a)]+6((>(x)}. (11) 

The corresponding gauge covariant derivatives, D, D 2 
and D A respectively are constructed by including link 
variables in the usual way. The chromoelectric field is 
discretised by a clover term with plaquettes in the three 
space-time planes only 

with 

Qi(x) = U l {x)U t {x + i)u}{x + i)ul{x) 

+ U t {x)u} (x - z + t) U} (x - 1) U t (x - 1) 

+ Uj(x - i)Ul(x - 1 - i)U % {x - i - i)U t (x - i) 

+ u}{x-i)u t {x- i) Ut{x + i-i)u}{x). 

(13) 

u s and Ut are the mean-link improvement parameters. 
u s is determined from the spatial plaquette and u t is set 
to unity. At the accuracy of the action constructed here 
no improvement is required. Finally, including the gauge 
fields and the mean-link improvement factors the lattice 
fermion matrix is given by, 
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M AmA ip{x) 



— \ [ a r ma t 



18s ra 2 a \ 1 

— + r + -pe.fi, J>(x) - — (r - l0 )U t (x)^{x + t) + (r + l0 )Ul{x - t)^{x - t) 
? 4 / 2u t 



i r i 2 12 

— y — (4s - -fi r Ji)Ui(x)ip(x + i) H (4s + -Hrli)UHx - i)ip{x - i) 

?<3 i L U s 3 u s 3 

- 4^)^^)^^ + t)if)(x + 2i) - ^(s + ^rli)U}{x - i)l//(a: - 2i>(a: - 2i) 1 . (14) 
u s 12 12 J J 



r 



At the tree-level, the fermion anisotropy £ g is given by the 
ratio of scales, £ = a s /a t . We call the action described 
here ARIA for Anisotropic, Rotated, Improved Action. 
It is classically improved to 0(a t , a^). 



III. HEAVY QUARKS WITH ARIA 

The precision calculation of the glueball spectrum on 
coarse lattices pj suggests that heavy hadronic quanti- 
ties would also benefit from the anisotropic formulation. 
The correlation functions for heavy-heavy and heavy- 
light mesons fall off rapidly with time and it can be 
difficult to isolate a convincing plateau over a reason- 
able number of timeslices. A lattice with fine tempo- 
ral direction in principle solves this problem by provid- 
ing a large number of timeslices over which the time- 
dependence can be resolved. Improved Wilson actions 
on anisotropic lattices have been used to study a range 
of heavy flavour physics including charmonium and bot- 
tomonium spectrosc opy MIH El- heavy-light and hy- 
brid spectra El El El 13 and also heavy-light semilcp- 
tonic decays \LT\ . 

In these calculations, currents are improved using rota- 
tions, which are applied identically in all four space-time 
directions and the Wilson parameter in the spatial direc- 
tion is usually chosen to be either r s — 1 /£ |3 or 
r s = 1 = r* |20|, |2lJ, |22j, |23j. However, it was pointed 
out in Ref. [5| that simulations with anisotropic Wilson- 
type actions may include 0(a s m q ) effects. Naively, errors 
of this form are unexpected but they arise in products 
of the Wilson and mass terms in the action. In par- 
ticular, the authors showed that the presence of these 
artefacts, which appear in radiative corrections, depends 
on the spatial Wilson parameter, r s . The 0{a s m q )- 
dependence potentially spoils the benefits of working on 
an anisotropic lattice, especially at large quark masses. 

In Ref. [|| a different approach was adopted. Since 
the unwanted C(a s mg)-dependent terms arise from the 
spatial Wilson term the authors propose an anisotropic 
D234-type action 0] may be more suitable. In this case 
a rotation term is applied in the temporal direction only, 
removing the temporal doublers. Spatial doublers are 
removed by adding an irrelevant, dimension-four term 
to the Dirac operator. The authors showed to one-loop 
order in perturbation theory, that this so-called "sD34" 



action does not suffer from 0{a s m q ) terms. Comparing 
the ARIA action proposed in Section ITT1 and the sD34 
action from Ref. Q we see that these are the same, up 
to 0(at) improvement. 

The D234 quark action on an anisotropic lattice @| is 
be written 

= a t a s 5^(aOM^(s), (15) 

X 

and writing M in the notation of Ref. [|| 
M = m + ^^7M^(l-^<A M )-ia t (^rA M 

+ Y c sw a ^ F ^ ] + Y ^ d ^ a l A l- ( i6 ) 

The sD34 action is a special case of this action in which 
the parameters have the following values 

(y 0l Vi) = (l ) v); (& ,&i) = (0,±); (d , d*) = (0, ±); 
(ro,n) = (r t ,0); (cg w ,c l sw ) = (0,c sw )- 

(17) 

and v = (1 + ir t a t TOo)- Substituting in Eq. Ijltifl gives 
M sD34 = to + Y f 1 _ fi°»^) + 70 V ° 

i \ / 

_££ ( rAo + clw^o) + \al Y vA l ( 18 ) 

i 

which is the action we use in our simulations, up to 0(a t ) 
improvement. Reexpressing the fermion matrix in our 
notation, 

Maria = frmo + Y^^ l ^ i \ l ~ ^ a ^ 4 ) + 70 V ° 

i 

-?l(rA - r Jla i0 F i0 )+ S a 3 s Y*l (™) 

i 

where s = 1/8 and \i r = (l + ^r t atTOo). 

A. Analytic results for ARIA 

In this section the energy-momentum behaviour of the 
ARIA action is calculated. We begin by presenting re- 
sults for general r and s. The free-quark dispersion rela- 
tion is obtained by solving det Maria = in momentum 



4 



space where Maria is the Fourier transform of A/aria 
in Eq. I|14|) . The energy-momentum relation is 



p = HrPi(l + -a 2 s p 2 ), 



(22) 



cosh(i?a t ) = 



ruj(p) 



r 2 - 1 



± 



1 2 
with pi — — sin(a s pi) and pi — — sm(a s pi/2). Expand- 
a. a.. 



ing the physical solution in powers of spatial momentum 

y/(r + uj{p)) 2 + (1 - r 2 )(l + aW) (2Q) yields 



r 2 - 1 



where to(p) and p are defined as 

uj(p) = a t p r m + a t 5 a 3 s pf 



E 2 (p)=M 2 + ^p 2 + 0(p% 

Mi 

(21) where M\ is the rest mass, given by 
I 



(23) 



Mi = — cosh" 
a t 



-l ( r 2 + fi r m ra t — \/l + 2^ r m ra t + m\[i 2 c 



(24) 



The kinetic mass, M2 is given by 



1 

~M~ 2 



Pr a t 



yi~+ 2fj, r m rat + 



miPv a t 



r 2 + Urm^rat - \/l + 2fi r m Q ra t + m\\i 2 a 2 t 



r 



(25) 



Eqs. Il24|l and 125|l indicate that at the tree-level M\ and 
M2 do not depend on 0(a s m q ) terms or on the ratio of 
scales, £. 

To compare these expressions with the results of other 
studies, the particular choice r = 1 was considered. In 
this case the lattice ghost (the unphysical solution of 
Eq. {2D}) disappears and the dispersion relation is given 
by 



4 sinh 



with 



Eat 
2 



1 



a 2 p 2 



u 2 (p) 



1 + w(p) 



Mi = — log(l + fj, r m Q a t ), 
a t 

1 [i r 



2M, 



m (2 + /J r moo t ) 



(26) 

(27) 
(28) 



where now fi r = (1 + \at mo). These expressions are 
consistent with those obtained in Ref. |(| for the sD34 
action and in Ref. 01 f° r the Fermilab action on an 
isotropic lattice. 

The free-quark dispersion relations for massless and 
massive quarks are shown in Figure ^ The anisotropy 
parameter, £ is six for both cases. In analogy to the 
traditional Wilson r-parameter, the parameter s in this 



action can in principle take any positive value. We chose 
.s = 1/8 by eye, demanding that the energy- momentum 
relations do not have negative slope for a s |p| < ir. Since s 
parameterises a term which removes the spatial doublers 
and is irrelevant in the continuum limit precise tuning is 
not required. 



IV. RESULTS 

In this exploratory study the temporal rotations have 
been omitted which leads to an (D(at) classical discretisa- 
tion error. However, since at is small in these simulations, 
a t ~ 0.04 fm, the effects should be under control at least 
when atm q < 1. Discarding temporal rotations means 
the action has no clover term and in addition we have 
set fi r = 1. It is planned to include correction terms to 
remove 0(at) errors in future work. 

The ratio of scales is changed in a simulation by quan- 
tum corrections. Therefore the action parameter must 
be adjusted so that the ratio of scales measured from a 
physical quantity is correct. In a quenched simulation 
the parameters £ g and £ g in the gauge and quark actions 
may be independently tuned to the target anisotropy, us- 
ing different physical probes. This is not the case for un- 
qucnchcd simulations where the anisotropy in the gauge 
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(1,0,0) axis 

1 - (1,1,0) axis 

(1,1,1) axis 




n -l2n -l3n 



aJPl 

FIG. 1: The dispersion relations given by Eq.H3with { = 6, 
r = 1 and s = 1/8. The top figure is the massless case while 
the bottom plot shows the massive case, with a t m q = 0.2. 



and quark actions must be tuned simultaneously [9|- 

For this study an ensemble of quenched gauge config- 
urations for which £ g had already been tuned was used. 
In this case the tuning criterion was that £ = 6 when 
measured from the static interquark potential in differ- 
ent directions on the lattice. The parameter £ q in the 
fermion action must now be tuned such that its value 
determined from the energy momentum dispersion rela- 
tion is six. At this point we introduce some terminology 
which makes clear the difference between £ g , which is a 
parameter in the action, and the slope of the dispersion 
relation which is a physical observable - usually called 
the speed of light, c. The target anisotropy is six. £ 9 
is tuned so that the speed of light (determined from the 
slope of the dispersion relation) is unity. 

The anisotropic action offers the possibility of preci- 
sion studies of a range of phcnomenologically interesting 
heavy quark quantities in the D, B, J/tp and T sectors. 
For this reason it is important to understand the depen- 
dence of £ q on the heavy quark mass used in simulations. 
In particular, a contribution of (D(a s m q ) to the renor- 
malised anisotropy would spoil this tuning for charm and 
bottom quark masses. The main result in this section is 
a study of the mass-dependence of the speed of light at 



fixed anisotropy. 

A. Simulation parameters 

The gauge action used in this simulation is a two- 
plaquette improved action designed for precision glue- 
ball simulations on anisotropic lattices. A description is 
given in Ref. |25j | . The construction of the fermion action 
is described in detail in Section [HI Details of the simu- 
lation and parameter values are summarised in Table [I] 
A broad range of quark masses was investigated, from 



# gauge configurations 


100 


Volume 


10 3 x 120 


a s 


0.21fm 


a s /r 


0.4332(11) 


£ = a s /a t 


6 


a t m q 


-0.04,0.1,0.2,0.3,0.4,0.5,1.0,1.5 



TABLE I: Details of the simulation. 



a t m q — —0.04 which is close to the strange quark on these 
lattices to heavy quarks with atm q = 1.0 and 1.5. Both 
degenerate and nondegenerate combinations are consid- 
ered. The nondegenerate combination is made with the 
lightest quark and each of the heavier quarks. Note that 
atm q = —0.04 corresponds to a positive quark mass since 
Wilson-type actions have an additive mass renormalisa- 
tion. We accumulated data at spatial momenta (0,0,0), 
(1,0,0), (1,1,0) and (1,1,1), in units of 2tt/ci s L, averaging 
over equivalent momenta. 

It is worth noting that all the gauge configurations and 
quark propagators used in this study were generated on 
Pentium IV workstations. Generating the lightest quark 
propagators (close to the strange quark mass) required 
approximately one week on a single processor. At this 
quark mass no exceptional configurations were seen. 

B. Effective masses 

The success of anisotropic lattice methods is predom- 
inantly due to the increased resolution in the temporal 
direction. The fineness of the lattice in this direction is 
particularly useful when determining heavy mass quan- 
tities whose signal to noise ratio decreases rapidly. The 
increase in resolution also leads to reduced statistical er- 
rors in effective masses since fits can be made to longer 
time ranges than is usually possible with an isotropic lat- 
tice. For the same reason, the fitted values tend to be 
less sensitive to fluctuations of one or two points in the 
chosen fit range. 

In this study the effective masses were determined us- 
ing single cosh fits with a % 2 minimization algorithm. 
The signal to noise ratio was enhanced by using four 
sources, distributed across the lattice at timeslices 0, 30, 
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60 and 90. The average of these results was used in the 
effective mass fits. The statistical errors shown are cal- 
culated from 1000 bootstrap samples in each fit. Fig- 
ure |21 shows four effective mass plots. The first plot is 
the pseudoscalar meson with degenerate quarks at the 
lightest mass for zero momentum and for momentum of 
(1, f , f ) in lattice units, 2 / n/a s L. The second plot is the 
analogous case for the degenerate combination of quarks 
with dfTtiq = f.0. In all cases a clear plateau, over a 
large number of timeslices is observed. The fits to ef- 
fective masses of the non-degenerate mesons are equally 
good and in all cases the fit range is ten or more times- 
lices with a x 2 P er degree of freedom (x 2 /Ndf) ~ 1. In 



0.6 



0.5 - "> 



a,m q = -0.04 
O n - (0,0,0) 
>n = (1,1,1) 




t/a, 



1.65 



O n p = (0,0,0) 
On =(1,1,1) 



1.45 



30 

t/a, 



FIG. 2: Pseudoscalar meson effective mass plots. The two 
plots indicate that very good fits can be made for a wide 
range of quark masses and momenta. The top plot shows the 
effective mass of the lightest meson made from a degenerate 



combination of quarks with a t m q 



-0.04 for zero momentum 



and for three units of momentum in lattice units. The second 
plot is the analogous case for atm q = 1.0. 



Figure |3 the equivalent results for vector mesons are pre- 
sented. Once again, the lightest and heaviest degenerate 
combinations of quark masses considered are shown and 
very good fits are possible in both cases. 
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FIG. 3: Vector meson effective mass plots. As in the pseu- 
doscalar case shown in Figure H good fits are achieved over 
a large number of timeslices for all the quark masses consid- 
ered in this study. The top plot shows the lightest degener- 
ate vector while the bottom plot shows the same result for 
atm q = 1.0. 



C. The renormalised anisotropy 

The slope of the energy momentum dispersion rela- 
tion is determined nonperturbatively and compared to 
the target anisotropy, £ = 6. We are interested in 
both the precision of the determination and the devia- 
tion of the speed of light from unity. The wide range 
of quark masses used in this simulation (fltm q = —0.04 
to atvn q — 1.5) allow us to study the mass-dependence 
of the renormalisation. We also examine the difference 
between the anisotropy determined from particles with 
degenerate and non-degenerate quark content. 

To begin, the dispersion relation was determined for a 
pseudoscalar meson made from the lightest quarks in this 
simulation, a t m q = —0.04 and with an input anisotropy, 
£ g = 6.0. The value of c determined from the disper- 
sion relation was used to determine the tuned value of 
the anisotropy, £ g = 6.17 and the calculation repeated. 
The resulting dispersion relation is shown in Figure ^ 
The subsequent value of c determined from this data is 
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1.02 ±0.01. This value of ^ g = 6.17 was then used in sim- 



0.05 - 



0.04 - 



0.01 " 




FIG. 4: The energy-momentum for the lightest degenerate 
meson in this simulation. The bare quark mass is atm q = 
-0.04. 

ulations for a range of quark masses, 0.1 < atm q < 1.5. 
A representative sample of the energy-momentum disper- 
sion relations for this range of quark masses and particles 
is shown in Figure |3J The plot shows very good linear 
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2S 
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(a s p) 2 



FIG. 5: The energy momentum dispersion relation for the 
pseudoscalar meson at all the masses simulated. The lightest 
point, atrriq = —0.04 is close to the strange quark while the 
heaviest mass is close to the bottom mass. The points have 
been shifted about their momentum value to make the plot 
easier to read. 

dispersive behaviour. This relativistic dispersion relation 
persists for both degenerate and non-degenerate quark 
combinations in pseudoscalar and vector particles at all 
masses. 

The mass-dependence of the speed of light is given by 
the difference in the slopes for the different masses. Fig- 
ure shows that this dependence is mild. The lightest 
mass, close to the strange quark mass, is the noisiest and 
the statistical errors increase with increasing momentum, 
as expected. It should be noted that the quark propaga- 



tors used in this study are generated with point sources 
and the use of smearing techniques is expected to improve 
the signal for this and lighter quark masses. In addition 
the advantages of stout link gauge backgrounds [2(| will 
be investigated in further studies. 

In Tables ITT1 and lllll we show the speed of light deter- 
mined from the slope of the dispersion relation for each 
mass in the simulation. The x 2 /Nd.f. for these fits also is 
shown. Results for both pseudoscalar and vector mesons 
are given and the ground state masses extracted in the 
fitting procedure described in Section IIV Bl are listed. 





Pseudoscalar 




Vector 




a t m q 


a t M PS 


c 


X 2 /N df 


a t M v 


c 


X 2 /N df 


-0.04 


0.1045^5 


1.02li 


6.3/2 


O.I6II2 


O.97I2 


0.66/2 


0.10 


0.38311* 


0.983l 7 


2.8/2 


0.39341* 


0.9821* 


2.1/2 


0.20 


O.54I8I4 


0.995l 7 


0.33/2 


0.54721* 


0.990lg 


2.1/2 


0.30 


0.68871* 


1.010l 7 


2.4/2 


0.69241* 


O.997I9 


4.5/2 


0.40 


0.8269+.4 


1.022+j 


0.65/2 


0.82941* 


l.Oll!;! 


2.3/2 


0.50 


0.956914 


1.035l| 


1.3/2 


0.95871* 


1.0251;? 


1.6/2 


1.00 


1.5086l| 


I.O69I5 


1.3/2 


1.50921^ 


1.072lg 


1.2/2 


1.50 


I.9428I3 


1.075±| 


0.081/2 


1.9431!! 


1.0721J 


0.058/2 



TABLE II: The ground state pseudoscalar and vector masses 
with degenerate quarks. The speed of light determined from 
the dispersion relation for each quark mass is shown with the 
associated \ 2 /Nd.f. ■ The errors in all cases are statistical only. 
The parameter, £ q is fixed in these simulations to 6.17, its 
value determined from the dispersion relation of the lightest 
degenerate pseudoscalar meson. 



The dependence of c on the quark mass in the simulation 
is shown in Figures and [7| The plots show the speed 
of light as a function of the meson mass in units of a t 
for both pseudoscalars and vectors. It is important to 
remember that the anisotropy was tuned only once at 
the lightest pseudoscalar particle. The plots show good 





Pseudoscalar 




Vector 




a t m q 


a t Mps 


c 


X 2 /N df 


a t Mv 


c 


X 2 /N df 


0.1 


0.26101s 


0.981} 


0.23/2 


0.2802l| 


0.981a 


0.19/2 


0.2 


0.3466lg 


1.01±* 


0.56/2 


0.360llg 


O.99I2 


0.64/2 


0.3 


0.4254l 7 


1.02+j 


2/2 


0.435llg 


l.OOlj 


0.45/2 


0.4 


0.4987l 7 


i.oilj 


1.5/2 


0.5056lg 


O.99I2 


1.4/2 


0.5 


0.5668lg 


1.02!^ 


1.7/2 


0.5720lg 


l.OOla 


1.6/2 


1.0 


0.852lli£ 


l.OOlj 


2.6/2 


0.854li 


1.02l^ 


0.62/2 


1.5 


1.074li 


1.02+J 


2.1/2 


1.075±i 




1.8/2 



TABLE III: The ground state masses of non-degenerate com- 
binations of quark masses. In each case the quark mass 
given is combined with the lightest mass in our simulations, 
atm q = —0.04. As in Table ITT1 the pseudoscalar and vector 
meson states are shown with the speed of light and the asso- 
ciated x 2 '' /N d f . Once again all errors are statistical only and 
£, = 6.17. 
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FIG. 6: The mass-dependence of the speed of light determined 
from the pseudoscalar (PS) dispersion relations for fixed £<j = 
6.17. The plot shows mesons with both degenerate and non- 
degenerate quark mass combinations plotted as a function of 
the degenerate meson mass, in units of the temporal lattice 
spacing. 
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FIG. 7: The mass-dependence of the speed of light determined 
from the vector (V) dispersion relations with £ q = 6.17. Both 
degenerate and non-degenerate mass combinations are shown, 
plotted as a function of the degenerate vector mass as in Fig- 
ure E| 



agreement between determinations of c from degenerate 
and non-degenerate particles up to atm q ~ 0.5, corre- 
sponding to a t Mps — 0.957(5) in Figure El The charm 
quark mass on this lattice is close to atm q = 0.2, imply- 
ing that charm physics is both computationally feasible 
and requires little parameter tuning at an anisotropy of 
six. Figures IH1 and also show some quark mass depen- 
dence for degenerate mesons with a t m q > 0.5. They also 
indicate that the agreement between the degenerate and 
non-degenerate meson physics decreases for atm q > 0.5. 
This is not unexpected since degenerate mesons with two 
heavy quarks (charmonium and bottomonium) have a 
small Bohr radius and can suffer large discretistation ef- 
fects. The non-degenerate mesons (D and B mesons) do 



not have such a problem. 

We have investigated this dependence by varying the 
parameter £ q in the quark action and repeating the 
simulations described above, for the heavy quark mass 
a t m q = 1.0. The dependence of the speed of light, 
determined from the dispersion relation, on the input 
anisotropy is shown in Figure [H] The value of c deter- 
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FIG. 8: The speed of light as measured from the dispersion 
relation as a function of the parameter, £ q in the action. The 
first plot shows the result for the pseudoscalar (PS) meson, 
for both degenerate and non-degenerate combinations. The 
second plot is the analogous result for the vector (V) mesons. 



In both cases the quark mass is fixed at atm q 



1.0. 



mined from the degenerate meson moves closer to its 
target value of unity and c determined from the non- 
degenerate physics moves away from this value. It is 
also interesting to note the agreement between determi- 
nations of c from pseudoscalar and vector particles. The 
tuning, described above at atm q = —0.04 was carried out 
for pseudoscalars and it is reassuring that although the 
vector particles have larger statistical errors they never- 
theless yield a consistent picture for the mass dependence 
of the speed of light. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have explored the viability of 
anisotropic actions for heavy quark physics. An action 
suitable for simulations at large anisotropies and which 
has no 0(a s m q ) errors is described. One of the main dis- 
advantages of using anisotropic actions is the extra pa- 
rameter tuning required to recover Lorentz invariance. In 
particular, if the ratio of scales £ is sensitive to the quark 
mass in the simulation then a parameter tuning may be 
required for each mass. We have measured the speed of 
light for a range of quark masses having fixed the ratio of 



9 



scales at the strange quark, atm q = —0.04. Only slight 
mass dependence (for the degenerate mesons) is found 
up, to a t m q — 0.5 which is heavier than the charm quark 
on these lattices. This implies that one measurement of 
the speed of light is all that is required for simulations 
over a large range of masses, at the percent-level of sim- 
ulation. The simulations were repeated for mesons with 
non-degenerate quarks, using a value of £ tuned from 
the degenerate meson spectrum. The results are in ex- 
cellent agreement up to a t m q = 0.5. Since the charm 
quark on this lattice is approximately a t m q = 0.2 this 
work indicates that both heavy-heavy (degenerate) and 
heavy-light (non-degenerate) charm physics can be eas- 
ily reached using an appropriately improved anisotropic 
action. 

The results also show that heavy-light as well as heavy- 
heavy physics can be reliably simulated after a single 
tuning of £ 9 . The determination of c can be inter- 
preted as a measure of the ratio M1/M2 in Eq. (|23|l . 
The agreement of Mi and M2 for both heavy-heavy 
and heavy-light systems can in turn be interpreted as 
an absence, in this quark action, of the anomaly first 
discussed in Ref. [22j|. This anomaly was explained in 



Ref. [28j where it was pointed out that for a sufficiently 
accurate lattice action (0(v A ) in NRQCD) the discrep- 
ancies in binding energies SB = B2 — B\ vanishes and 
/ = (2SB Qq - (SBqq + 5B qq ))/2M 2Qq = as expected. 
The action described in this study has this property. 

This study has been carried out in the quenched ap- 
proximation which is a useful laboratory in which to 
study mass-dependent and tuning issues at relatively low 
computational cost. We are currently developing algo- 
rithms for dynamical simulations with anisotropic lat- 
tices which we plan to use in a study of heavy-flavour 
physics. 
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